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ON  STRONG  SLOT  INJECTION  INTO  A  SUBSONIC  LAMINAR  BOUNDARY  LAYER 

M.  Napolitano  and  R.E.  Messick 

ABSTRACT 

This  paper  is  concerned  with  the  problem  of  strong  slot  injection  into 
a  subsonic  laminar  boundary  layer  at  asymptotically  high  Reynolds  number. 

The  problem  is  formulated  and  the  governing  equations  are  presented  within 
the  context  of  triple  deck  theory.  The  linear  problem,  valid  when  the 
injection  velocity  is  small,  is  solved  analytically  by  Fourier  Transform 
techniques.  Graphical  results  are  given  for  a  wide  range  of  slot  lengths. 

A  numerical  technique  for  the  nonlinear  equations  is  presented  and  used  to 
obtain  solutions  for  various  injection  velocities.  Separation  is  found  to 
first  occur  downstream  of  the  slot,  where  a  recirculating  flow  bubble  is  formed. 
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I.  INTRODUCTION 


Injection  of  a  secondary  fluid  into  a  well  developed  boundary  layer  is 
a  widely  used  technique  in  the  aeronautical  industry.  It  is  used,  for  example, 
for  cooling  gas  turbine  blades  and  for  controlling  transition  and/or  separation 
(blow  off)  of  the  boundary  layer  over  airplane  control  surfaces.  It  is  therefore 
of  great  practical  significance,  for  the  aerodynamicist  and  the  thermal  analyst 
as  well,  to  be  able  to  analytically  predict  the  effect  of  injection  on  the  flow 
pattern  and  temperature  distribution  around  a  turbine  blade  or  any  other  piece 
of  machinery.  Our  interest  is  presently  limited  to  the  aerodynamic  aspect  of 
the  problem.  In  this  respect,  most  of  the  earlv,  classical  work  on  injection 
is  concerned  with  the  case  of  weak  blowing  into  an  incompressible,  laminar  boundary 
layer;  that  is,  the  ratio  between  injection  and  free  stream  velocities  is  assumed 
to  be  inversely  proportional  to  the  square  root  of  the  characteristic  Reynolds 
number.  Re.  For  this  case,  classical  boundary  layer  theory  is  valid  and  self 
similar  solutions  are  possible  [1],  The  most  important  conclusion  of  such 
studies  is  that  weak  injection  has  a  much  greater  effect  on  the  boundary  layer 
properties  (wall  shear,  displacement  thickness  eta)  than  one  might  expect,  and 
such  Injection  can  even  blow  the  layer  off  the  plate,  i.e.  induce  separation. 
Subsequently,  the  interest  of  most  researchers  has  been  attracted  by  the 
supersonic  mainstream  case.  See  Smith  and  Stewartson  [4]  for  a  critical  review 
of  their  efforts.  In  particular,  after  the  formalization  of  the  asymptotic  triple 
deck  theory  of  Stewartson  et  al.  [2,3]  a  rigorous  definition  for  moderate  and 
strong  blowing  has  become  available,  together  with  the  appropriate  governing 
equations.  Smith  and  Stewartson  [4]  consider  the  problem  of  strong  slot  injection 
into  a  supersonic  laminar  boundary  layer;  for  this,  they  provide  analytical 
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solutions  to  the  linear  problem,  valid  when  the  injection  velocity  is  small 
in  the  triple  deck  scaling,  and  a  numerical  technique  for  the  full  nonlinear 
problem.  Such  a  numerical  technique,  though,  developed  a  "pressure  sensitivity" 
or  instability,  due  to  ineffective  imposition  of  the  downstream  boundary 
condition,  and  could  not  account  for  separation  ahead  of  the  slot.  Napolitano 
[6]  has  recently  overcome  these  two  difficulties  and  provided  more  complete 
results,  in  which  separation  is  indeed  found  to  first  occur  ahead  of  the  slot. 
Finally,  Smith  and  Stewartson  [5]  also  considered  the  problem  of  strong  (and 
massive)  plate  injection  into  a  separated  supersonic  laminar  boundary  layer. 

Sychev  [7]  has  recently  shown  that  triple  deck  theory  is  valid  for  subsonic 
flow  conditions  as  well.  The  problem  of  strong  slot  injection  into  a  subsonic 
laminar  boundary  layer,  has  thus  become  feasible,  and  is  the  subject  of  this 
paper.  In  section  II  the  problem  is  defined  and  the  triple  deck  governing  equations 
presented.  In  section  III,  the  case  in  which  the  injection  velocity  is  small, 
is  considered:  the  governing  equations  are  linearised  and  solved  analytically 
by  Fourier  transform  Techniques.  In  section  IV  a  numerical  technique  for  the 
full  nonlinear  problem,  derived  from  that  of  Napolitano  et  al  [8],  is  presented. 
Finally,  in  section  V  results  are  provided  for  the  linear  and  nonlinear  problems. 
Linear  pressure  and  wall  shear  solutions  are  presented  for  a  wide  range  of  slot 
lengths.  Nonlinear  results  are  instead  provided  for  one  slot  length  only  but  for 
various  values  of  the  injection  velocity. 


II.  Problem  Definition  and  Governing  Equations 


A  Newtonian  fluid,  of  kinematic  viscosity  v*,  flows  past  a  semiinfinite 

flat  plate  with  subsonic  free  stream  velocity,  U*,  parallel  to  the  plate.  At 

00 

a  distance  x*  from  the  leading  edge  of  the  plate,  fluid  is  injected  perpendicular 

to  the  wall,  with  velocity  v^*,  thru  a  slot  of  length  x*  -  x*.  The 

characteristic  Reynolds  number.  Re  =  IT*  x*/v*,  is  asymptotically  large  and  the 

—3/8  —3/8 

injection  velocity  and  slot  lengths  are  of  order  Re  U*  and  Re  x*, 

oo  O 

respectively.  For  this  flow  situation,  (see  Figure  1)  three  region  or  decks 

develop  around  the  injection  region,  the  lower,  middle  and  outer  decks,  all 

3 

having  the  characteristic  length  of  the  interaction  region,  of  order  e  x* 

o 

(e  =  Re  and  characteristic  height  of  order  e^x*  ,  x*  and  x* 

o  o  o 

respectively.  The  middle  deck  is  a  displaced  31asius  boundary  layer,  having 
the  passive  role  of  communicating  the  pressure  interaction  between  inner  and 
outer  decks.  This  last  one  is  a  potential  flow  past  a  "thin  airfoil" 
produced  by  the  viscous  effects  in  the  lower  deck  underneath.  The  inner  (lower) 
deck,  finally,  is  a  viscous  region  which  arises  from  the  need  to  satisfy 
the  wall  boundary  conditions  on  the  plate.  After  appropriate  scaling  of  the 
independent  and  dependent  variables  and  proper  matching  between  the  neighboring 
decks,  the  Navier  Stokes  equations  reduce  to  the  triple  deck  ones.  These 
are  basically  the  lower  deck  boundary  layer  like  equations  coupled  to  the  outer 
deck  (linear  airfoil)  pressure  interaction  law.  The  continuity  and  momentum 
equations  are  given  as: 


u  +v„*0;  uu  +  vu  “ 
x  y  x  y 


.  Me. 

dx 


+  u 


yy 


(2.1,a,b) 
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with  boundary  conditions: 


u(x,0)  =  0 


(2.2) 


v(x,0) 


for  0  <  x  <  x 
—  s 

elsewhere 


(2.3) 


u(x,y  -*■«)-*■  y  +  6, 


(2.4) 


u(x  -*■  ,y)  -*■  y  . 


(2.5) 


The  pressure  displacement  interaction  law  is  finally  given  by 


6 

xx 


p?d£ 

* 


(2.6) 


where  the  integral  denotes  the  Cauchy  principal  value.  Note  that  5  is  the 
total  displacement,  due  to  the  viscous  effects  and  to  the  injected  fluid  as 
well.  Also  note  that  6  and  p  are  independent  of  the  normal  coordinate  y. 
For  these  two  quantities  the  x  subscript  indicates  a  total  (rather  than  a 
partial)  derivative. 
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III.  Solution  to  the  Linear  Problem 


The  problem  described  by  equations  (2.1)  thru  (2.6)  is  nonlinear,  because 
of  the  convective  terms  u  ux  and  v  u^  in  the  momentum  equation  (2.1b)  and, 
at  present,  can  be  solved  only  numerically.  The  case  in  which  the  injection 
velocity  Vw  is  very  small  (Vw<<1),  however,  can  be  studied  analytically.  This 
analytical  solution  is  very  important  since  it  provides  interesting  information 
on  the  asymptotic  decay  of  the  solution  for  x  -*■  +  °°  and  is  the  only  tool 
for  assessing  accuracy  and  reliability  of  the  numerical  results.  Moreover,  its 
usefulness  can  be  reasonably  expected  to  be  even  greater,  because,  for  the  super¬ 
sonic  mainstream  case,  it  has  been  found  that  the  linear  solution  is  relatively 
close  to  the  full  nonlinear  one,  for  Vw  as  large  as  1  [6]  .  A  broad  range  of 
validity  for  the  linear  solution  has  also  been  found  for  both  supersonic  and 
subsonic  flow  past  a  hump  on  a  flat  plate  [8,9].  For  Vw  equal  to  zero,  the  flow 
field  is  one  characterised  by  uniform  shear  (u=y)  with  v,  p  and  5  identically 
equal  to  zero.  For  Vw<<1,  thus,  the  solution  is  a  small  perturbation  of  such  a 
uniform  shear  flow,  that  is 


u(x,y)  -  y  +  V u  (x,y)  +  ...  (3.1) 

w 

v(x,y)  -  Vwv(x,y)  +  ...  (3.2) 

<5(x,y)  »  V  o(x,y)  +  ...  (3.3) 

w 

p(x,y)  -  V  p(x,y)  +  ...  (3.4) 

w 
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I 

I 

In  order  to  obtain  the  governing  equations  for  the  linear  problem,  equations  | 

(3.1  -  3.4)  are  substituted  into  equations  (2.1  -  2.6)  and  only  the  terms 

containing  V  linearly  are  retained,  to  give: 
w 

u  +  v  *  0;  y  u  +  v=*  -  +  u  (3.5  a,b) 

x  y  x  dxyy 


u(x,0) 

v(x,0) 

u  (x,y 


=  0 


oo)  * 


0  <  x  <  x 
-  s 

elsewhere 

& 


(3.6) 

(3.7) 

(3.8) 


u  (x  -*■  -<=°,y)  =  0 


(3.9) 


x-e 


(3.10) 


Linearization  makes  it  possible  to  construct  the  solution  for  the  case 
of  a  finite  slot  from  the  solution  for  the  infinite  slot  (xg“»)  by  translation 
and  superposition.  For  clarity  and  convenience  the  superior  bars  will  be  omitted 
from  the  dependent  variables  for  the  case  of  the  infinite  slot. 

The  infinite  slot  problem  may  thus  be  formulated  as  follows: 


u  +  v 
x  y 


y  u  +  v  - 


-  p  + 
rx 


yy 


(3.11  a,b) 


i 

i 

i 
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for  (x,y)  e{ (-",“)X[0,«) } ,  with  the  subsidiary  conditions: 


u(x.O)  -  0 

(3.12) 

v(x,0)  «  h(x) 

(3.13) 

u  (-®,y)  ■  0 

(3.14) 

u(x,oo)  *  6 (x) 

(3.15) 

p(+o°)  «  0 

(3.16) 

i  °VP,(C)<3€ 

6  (x)  -  -  -  1  > . .  - 

XX  IT  f  X-C 

(3.17) 

where  h(x)  denotes  the  Heaviside  step-function  and  the  integral  the  Cauchy 
principal  value. 

The  system  (3.11)  may  be  transformed  to  a  single  equation  by 
differentiating  (3.11  b)  with  respect  to  y  and  using  (3.11  a)  to  eliminate 
v  .  The  result  is 

y 

V  "  y  w*  *  0  (3-18) 

where  w  *  u^.  Let  capital  letters  denote  the  Fourier  transforms  of  the 
corresponding  (lower-case)  dependent  variables  e.g. 

00 

W  -£  w  e*lux  dx  (3.19) 
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The  Fourier  transform  of  (3.18)  using  (3.14)  then  provides  the  Airy 
differential  equation 


W  -  itoy  W  »  0 

yy 


(3.20) 


where,  again  for  clarity  and  convenience,  the  product  ito  will  be  understood  to 
stand  for  e  +  ito,  with  e  >  0  an  arbitrarily  small  constant  to  assure  convergence 
of  the  transform  at  the  upper  limit.  The  limit  as  e  +  0  is  to  be  used  in  all  results. 
The  only  non- integral  powers  of  ito  to  be  encountered  here  involve  the  cube  root.  A 
natural  branch  for  these,  leaving  the  real  to  axis  free  is 


Lin 

e+0 


1_ 

[e  +  ito]3 


(3.21) 


for  — J  which  defines  the  cut  indicated  in  the  figure. 


The  general  solution  of  (3.20)  is  usually  written  as 


_1_  l_ 

W  -  CQ  Ai  ([ito]3  y)  +  CL  Bi  ([id)]3  y) 


(3.22) 


where  Cq,  are  constants  and  Ai,  Bi  are  the  Airy  functions  of  the  first  and 
second  kind,  respectively.  However,  it  now  follows  that  W  is  exponentially 
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unbounded  as  y  t  ®,  violating  conditions  (3.14)  and  (3.15),  unless  ■  0. 

With  *  0,  condition  (3.12)  implies  that 

y  — 

U  *  C  f  Ai  ([ia)]3OdC  (3.23) 

0  J0 

The  Fourier  transform  of  (3.15)  then  requires  that 

CO  1 

A  -  CQ  j£  Ai  ([ia)]3OdC  (3.24) 

However,  for  precisely  the  branch  defined  in  (3.21),  the  integral  in  (3.24) 

evaluates  to  - j-  •  Thus, 

3[iu]7 

1_ 

CQ  -  3[ioi]3’  A  (3.25) 

Now  A  is  related  to  P  by  means  of  the  transform  of  (3.17),  using  (3.16).  i.e. 


ituA  *  i  P  sgn  w 


where  sgn  w 


ai  >  0 

u  <  0 


(3.26) 


Condition  (3.13)  remains  to  be  satisfied.  However,  by  continuity  it  can  be 
used  with  (3.11  b)  to  give: 
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(3.27) 


u  *  p  +  h(x)  for  y  +  0 
yy  x 


The  transform  of  this  (with  our  limiting  convention)  is 


U  ■  iuP  +  7—  for  y  +  0 

yy  i(o 


(3.28) 


Substituting  from  (3.28)  for  U,  we  have 


CQ[iio] 3  Ai'(O)  =  ia.P  +  ^ 


(3.29) 


Equations  (3.25),  (3.26)  and  (3.29)  maybe  solved  to  give 


2_ 

P  .  -1M  "3 

D 

4 

r  «  -3i(in>)  3sgn  no 
0  D 

5_ 

-Kirn)  3  sgn  a 


(3.30a) 


(3.30b) 


(3.30c) 


4_  4_  .3. 

where  D  *  (iw)3  +  03  i  sgn  u  and  0  ■  [-3Ai(0)]U  ■  0.8272  ...  is  Lighthill's 

constant. 

U  is  now  known  from  (3.30b)  and  (3.23),  so  its  inverse  transform  u,  together 
with  p  from  (3.30a),  maybe  used  with  (3.11)  to  find  v  and  check  the  results. 

The  inverse  transforms  can  be  put  into  a  suitable  real  form  after  some 
cumbersome  but  straight-forward  manipulations.  These  manipulations  involve 
altering  the  paths  of  integration  to  radial  lines  at  angles  of  <f>  ■  +  ^  in  the 
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co-plane,  depending  upon  whether  or  not  +  x  >  0. 

The  result  for  the  pressure  then  takes  the  compact  form: 


p(x) 


00 


P^-1/3  pSi(x)+l 


(3.31) 


It  now  follows  by  superposition  that  the  pressure  for  the  finite  slot 
Is: 


p(x)  -  p (x)  -  p(x-xg)  (3.32) 

It  is  obvious  that  p(x)  is  continuous  everywhere  except  possibly  at  the  slot 
ends  i.e.  at  x*0,  xg.  Let  the  jump  in  a  function  f(x)  be  defined  as 


(f(x)}  -  f(x+)  -  f(x-) 


(3.33) 


then 


(p(0)}  -  -{p(x)}  -  (p(0)} 


3 

OO  00 

if  dP  0  f 

2ir0 

'—'0  p8-^  p*+l  j  0 

CD 

i 

1 _ 

jLr 

"/2ttcos  F  „  /T'fcos  8 

1  ■  ° 

2ir6  L 

2  4 

•  — 

J 

(3.34) 
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Also, 


(Px(0)}  *  "{Px(xs)}  “ 


This  shows  that  p  Is  (i.e.  may  be  defined  to  be!)  continuous  on  (-“>,«>),  but 

that  Its  derivative,  the  pressure  gradient,  displays  unit  jump  discoutlnuities  at 

the  ends  of  the  slot  in  opposite  directions. 

Considering  (3.31)  and  (3.32),  it  is  seen  that  outside  of  the  slot  the 

8  r-  4 

denominators  of  both  integrands  in  (3.32)  are  the  6ame.  (Note  that  p  -  /3p  +  1>0 

for  all  p  so  p(x)  >  0  for  all  x.)  However,  for  x  <  0  the  exponential  function 

in  the  integrand  of  p(x)  exceeds  that  in  p(x-x  ).  The  situation  is  reversed  for 

s 

X  >  X  . 

s 

It  follows  then,  because  of  the  difference  in  p's  in  (3.32),  that  p(x) 

is  positive  to  the  left  of  the  slot  (x<0)  and  negative  to  the  right  of  the  slot 

(x>x8).  Also,  differentiation  of  p(x)  for  x  <  0  does  not  alter  this  sign  relationship 

while  for  x  >  x  it  reverses  it.  This  means  that  in  both  cases  tT(x)  is  a  monotone 
8 

increasing  function  of  x.  Now  p’(j^)  ■  0,  so  that  n(x)  starts  from  zero  at  x  ■  -» 
and  increases  positively  to  the  maximum  of  p(0)  >  0  at  x  *  0.  To  the  right  of  the 
slot  p(x)  starts  at  the  negative  value  p"(x  )  for  x  ■  x  and  increases  to  zero 
as  x  +  ®. 
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Within  the  slot,  differentiation  of  p(x)  results  in  the  derivative  of 

p(x)  being  negative  while  the  derivative  of  p(x-x  )  remains  positive.  This  means 

s 

that  the  derivative  of  p(x)  within  the  slot  is  negative  and  that  p(x)  is  monotone 


decreasing  from  p(0)  > 0  to  the  minimum  of  p(x  )  <  0  at  the  r.h.  end  of  the  slot. 

s 

Finally,  it  is  clear  that  for  x  <_  0,  p(x)  is  a  monotone  increasing  function 

or  x  .  The  maximum  p"(0)  increases  from  0  for  x  *  0  to  the  maximum 
s  s 

3  ^2*"  it  _ . 

P(0)  *  cos  -jj-  ,  (see  (3.34)),  for  xg+  «  .  Conversely,  for  x  xg,  p(x)  is  a 

monotone  decreasing  function  of  x  .  The  minimum  p(x  )  decreases  from  0  for  x  *  0 

s  s  s 

2/2  TT 

to  the  minimum  -p(0)  ■  -  7——  cos  for  x  t  “  . 

49  0  s 

The  skin  friction  coefficient  T(x)  *  u^(x,0)  is  also  of  primary  importance. 
For  the  infinite  slot  t(x)  »  u^(x,0)  can  be  written  in  the  compact  form: 


2  3 
2ir  u 


f  w(9x)T 

L  r(l) 


«>  2  -^xlo^”^  4  ^(x) 

+  f  P  e  B|X>P  (T  p  -1) 

Jq  p®  -  i^3p^h(x)  +  1 


h(x)  -» 
ii - (3.: 

+ 1  J 


So  that,  as  before,  for  the  finite  slot; 


T (x)  *  t(x)  -  T(X-X  ) 


(3.37) 


Here  too  t (x)  is  obviously  continuous  except  possibly  at  the  slot  ends. 
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{t(0)}  -  -{t(x  )}  **{x(0)} 
s 


9/3 

2  3 
2ir  6 


f”ir/2  sin 
[-4-^ 


tt/2  sin 


] 


(3.38) 


Also, 


T 

X 


(x) 


9/3 

,  2.2 
2ir  e 


1  h(x)  « 

(-ao  r 

J0 


-e!x|p3 

[p  -/3ph(x)]e  dp 
p8-/?p^h(x)  +  1 


(3.39) 


Therefore: 


{T  (0)}  -  -{x  (x  )}  *{x  (0) } 

X  X  s  X 


9/3 

2  2 
2ir  6 


(p5~/3p)dp 

p8-/3p4+1 


9/3  f  /2*  /5w  “1  „ 

“  ■  ~J"° 


(3.40) 


Again, 


T  (x) 
XX 


9/3 

2ir26 


1  h(x) 
<2> 


1-  r.(to 

»  — 0  ] x]  p'5  .  — . 

r  e  l|K  dp  1 

L36^(x)T 

Jq  p8-/3ph(x)  +1  J 

(3.41) 
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So  that 


t  GO 

xx 


9/3  r~Ah(x)  r(3)  /Tttcos  8 


.  -*n  rci 

,.K  L(2 


1  1 
36T(x)y 


as  |x|  \  0  (3.42) 


In  summary  then  t(x)  and  its  derivative  are  continuous  on  but 


t  (x)  -  0(1/ [  X t 1//3)  or  0(1/  |x-x  t1^3)  as  x  tends  to  0  or  x  ,  respectively. 


T 


IV.  NUMERICAL  TECHNIQUE 


The  numerical  technique  of  Napolitano  et.  al.  [8,9],  when  applied  to  the 
triple  deck  equations  for  subsonic  flow  past  a  parabolic  hump,  successfully 
accounted  for  the  discontinuous  pressure  gradients  at  the  leading  and  trailing 
edges  of  the  hump.  Therefore,  it  could  be  reasonably  expected  to  be  applicable 
to  the  present  strong  slot  injection  problem  by  changing  the  wall  boundary 
condition  for  the  normal  velocity  component.  However,  another  modification  was 
found  necessary,  i.e.  the  downstream  boundary  condition  for  the  displacement 
thickness  6  had  to  be  applied  as  a  derivative  one,  ^  ■  0,  in  order  to  obtain 
correct  results.  This  was  needed  because,  due  to  the  injected  fluid,  6  decays 
to  zero  very  slowly  as  x  +  When  5  was  set  equal  to  zero  at  a  finite  downstream 
location, very  inaccurate  results  were  produced,  probably  due  to  an  artlfical 
blockage  of  the  flow.  The  derivative  condition  was  instead  quite  satisfactory 
and  was  used  throughout  this  study.  In  the  results  section,  the  wall  shear 
distribution  obtained  by  this  technique  for  Vw  -  0.1  is  shown  to  be  in  good  agreement 
with  the  linear  solution.  However j  for  larger  values  of  the  injection  velocity  Vw>  the 
slower  convergence  rate  of  the  numerical  technique,  together  with  the  necessity  of 
imposing  the  upstream  and  downstream  boundary  conditions  farther  away  from  the  slot, 
caased  the  computational  time  to  become  prohibitive,  An  effort  was  made,  therefore  to 
remove  these  two  major  limitations.  A  stretching  transformation  of  the  longitudinal 
coordinate  x  was  used  (as  in  Ref.  6),  to  reduce  the  number  of  computational  grid 
points  without  affecting  the  accuracy.  Also,  a  direct  (non-iterative)  solution 
of  the  second  sweep  equation  was  employed  to  enhance  the  convergence  rate.  The 
numerical  technique  as  modified  for  the  present  problem  is  outlined  below,  details 
being  provided  only  for  the  new  features. 


A  relaxation-like  time  derivative  of  the  displacement  thickness  5  was 
added  to  the  right  hand  side  of  the  momentum  equation  (2.1b)  so  that 


uu  +  vu 
x  y 


dpx  ,36 
dx  yy  3t 


(4.1) 


Then,  a  two-sweep  Alternating  Direction  Implicit  (ADI)  technique,  which  is 

accurate  to  second  order  in  time,  was  used;  the  computation  proceeded  from  time 

n  ^  n  At  *  n+1  *  At 

t  to  time  t  ■  t  +  ~  and  from  t  to  t"  =  t  +  as  follows: 


(uu  +  vu  -  u  )  *  -  (^)  +  X(5*  -  6n) 

x  y  yy  dx 


(4.2) 


(uu  +  vu  -  u  )  =  -  (^2.)  +  X(6n+^  -  6  ) 

x  y  yy  dx 


(4.3) 


with  X  »  . 


(4.4) 


In  practice,  at  the  *  time  level  (^£)  and  6  are  known  at  all  the  x^  locations 

and  equation  (4.2)  is  solved,  coupled  with  the  continuity  equation  (2.1a)f by 

*  *  * 

finite  difference  methods.  These  provide  the  u  ,  v  and  5  arrays.  This  step 
is  identical  to  that  used  by  Napolitano  [6]  to  solve  the  supersonic  mainstream 
strong  injection  problem  and  those  details  will  not  be  repeated  here.  A  second 
sweep  equation  for  the  ADI  procedure  is  needed  to  obtain  the  new  time  level 
values  of  (^)  and  6  .  The  whole  process  can  then  be  repeated  until  a 

suitable  convergence  criterion  is  met.  Such  an  equation  is  obtained  from 
equations  (4.2),  (4.3)  and  the  interaction  law,  equation  (2.6)  as  follows:  since  ^ 
and  6  do  not  depend  on  the  normal  coordinate  y,  equations  (4.2)  and  (4.3)  were 
combined  so  as  to  eliminate  all  of  the  velocity  terms. 
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(4.5) 


-  <&>” +  i<5"+1  - 2s* +  s!,) 


Then,  this  expression  for  (^) 


n+1 


is  introduced  into  the  interaction  law 
equation  (26)  which  is  also  evaluated  at  the  n+1  time  level.  This  gives 

n 


A 

dx2 


M  ,  -  <£>  * 

•  -1  i  — ■ — =— 

TT  *  X-£ 


-  25  +5  ) 


dg 


(4.6) 


n+1  do 

Next,  equation  (4.6)  is  solved  for  6  and  the  pressure  gradient  (^p) 


n+1 


is  evaluated  explicitely  by  means  of  equation  (4.5),  which  completes  the  two 

sweep  ADI  process.  Note  that  equation  (4.6)  is  a  second  order  differential 

equation  and  can  satisfy  both  the  upstream  and  downstream  boundary  conditions 

on  6.  Together,  these  constitute  a  properly  formulated  two-point  boundary 

n+1 

value  problem.  In  order  to  obtain  the  6  array,  by  numerically  solving 
equation  (4.6)jthe  integration  domain  is  limited  to  a  finite  range  x^,  x^.  where 
x^  and  Xj  are  locations  sufficiently  far  upstream  and  downstream  of  the  slot, 
respectively.  Equation  (4.6)  is  then  modified  to 


,2.  "-1  *I.«1 

dx  x,  * 


-1  (dp/dg)n  +  X(5n  -  26*) 

IT  '  X-C 

*1 


dC  (4.7) 


In  References  (8,9 )it  was  shown  that  the  error  in  the  Cauchy  integral  due 
to  the  omitted  integration  (-«,x^)  and  (x^,*)  was  negligible  for  a  proper  choice 
of  x^  and  x^ .  In  the  present  study,  it  was  also  verified  that  moving  x^  and  x^ 
farther  from  the  slot  did  not  produce  any  significant  change  in  the  results.  In 
equation  (4.7),  the  Cauchy  integral  has  been  split  in  order  to  place  the  unknown 
terms  on  the  left-hand  side  of  the  equation.  The  Cauchy  Integrals  here  were 
reduced  to  finite  difference  form  by  the  methods  given  in  Reference  9,  i.e., 
the  6  and  (dp/dg)n  functions  were  taken  as  constant  over  each  mesh  of  the 
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computational  domain  and  removed,  so  that  the  singular  terms,  — -  ,  could  be 

integrated  analytically  from  each  grid  point  to  the  following  one.  This  provided 

d  n 

the  proper  coefficient  for  the  previously  extracted  djand  (-77)  values. 

d26  0+1  i 

Also,  ( — x)  was  expressed  in  finite  difference  form  using  the  variable  grid 

dx 

second  order  accurate  representation  of  Blottner  [10]  i.e. 


'i-1 


(xi+rxi-i) 


Xi"Xi-l 


-  M- 


1  Wi 


Xi+l“Xi 


)  + 


i+1 


Xi+1-Xi 


(4.8) 


As  described  above,  equation  (4.7)  is  reduced  to  a  system  of  linear  algebraic 
equations  of  the  form: 


(4.9) 


whose  solution  finally  provides  the  sought  after  6^  array.  The  main  difference 
between  this  technique  and  that  of  References  (8,9)  lies  in  the  solution  procedure 
for  equation  (4.9).  This  will  now  be  discussed. 

In  References  (8,9)  an  iterative  procedure  was  used  to  solve  the  system  (4.8) 
at  every  n+1  time  level.  Equation  (4.9)  was  first  rewritten  as 


[Tij  +  (Aij  "  Tij)]  6j+1  "  >  (4.10) 

where  T^  is  the  tridiagonal  matrix  containing  all  the  A^  elements  belonging 
to  its  main, upper  and  lower  diagonals.  The  (A^  -  T^)  contributions  were 

then  evaluated  at  a  previous  iteration  level  and  combined  with  the  known  right- 
hand  side  vector  to  give 

Tu  *r  ■  *i  •  (a«  *  v  c  (‘-u> 

This  was  solved  repeatedly  by  the  very  fast  Thomas  Algorithm  until  a  convergence 
criterion  for  and  was  satisfied. 

pj  * 
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for  all  at  the  beginning  of  the  computer  program.  The  inversion  is  obtained  by 
means  of  the  Fortran  double  precision  version  of  the  IMSL  LIB1-0006  LINV2F 
subroutine  designed  for  ill-conditioned  matrices  such  as  the  present  one. 

It  uses  the  Gauss  elimination  procedure  with  iterative  improvements  until  the 
inverse  is  correct  to  a  prescribed  number  of  significant  figures  (chosen 

here  to  be  eight) .  By  using  this  inverse,  the  solution  of  the  second  sweep 
equations  involved  only  a  very  fast  direct  matrix  multiplication,  i.e. 


(4.12) 


In  practice,  the  matrix  was  reduced  to  a  1-2  by  1-2  square  matrix,  I  being 

the  total  number  of  grid  points  in  order  to  properly  satisfy  the  boundary  con¬ 
ditions  on  6n+1.  These  conditions  were  taken  to  be  {(x^)  ■  0  and  (Xj)  ■  0, 
respectively,  which  in  finite  difference  form  became 


&1  -  0  ;  «x  -  «x-1  (4 ,13a, b) 

The  I  by  i  to  1-2  by  1-2  matrix  reduction,  needed  to  enforce  equations  (4.13a,b) 
is  very  straightforward  and  is  omitted  here.  Note  that  in  equation  (4.12)  the 
indexes  i  and  j  vary  between  2  and  1-1.  In  the  computer  program  they  were  all 
reduced  by  one  to  accommodate  the  Fortran  language  implementation. 

In  Reference  (9)  it  was  shown  that  for  a  constant  step  size  all  the  logar¬ 
ithms  needed  for  evaluating  the  Cauchy  Integrals  in  equation  (4.6)  need  be 

evaluated  only  once  and  stored  in  a  vector  of  dimension  21 »  In  this  way  the 

> 

computational  effort  was  reduced  drastically  with  respect  to  evaluating  all  the 
logarithms  at  every  second  sweep  of  the  solution  process.  In  the  present  study 
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due  to  the  variable  grid  distribution,  a  1  by  I  matrix  was  necessary  to  store 
all  the  above  mentioned  logarithms. 


In  the  results  section  it  can  be  seen  that  the  improved  algorithm  was 

able  to  compute  the  solution  for  V  *  1.5  in  a  reasonable  CPU  time.  The 

w 

technique  of  References  (8,9)  in  this  case,  would  have  required  a  significantly 
higher  computer  storage  space  and  a  CPU  time  one  order  of  magnitude  higher. 
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V.  RESULTS 


a)  Linear  Equations  Results 

A  small  computer  program  was  used  to  evaluate  the  linear  theory  pressure 

and  wall  shear  given  by  equations  (3.31,  3.32)  and  (3.36,  3.37).  All  the 

integrals  were  evaluated  by  the  simple  trapezoidal  rule  with  the  range  of 

integration  for  p  limited  to  0  <  p  <  p  «  100.  It  can  be  seen  that  this 

—  —  max 

-A 

approximation  leads  to  an  error  of  order  pmny  or  less,  which  is  obviously 
negligible.  However,  the  integrand  of  the  integral. 


*^3/2  p°-p 


6  2 


o  p 


8 


-  Jl  p*+l 


is  of  order  — ^  as  p+«. 


dp  (5.1) 

This  appears  in  the  expression  for  the  wall  shear 


tw»  equation  (3.36), when  either  x  or  x-xg  are  equal  to  zero.  It  Is  easily  seen 


that. 


max 


J 

o 


Jin 


6 

p  - 


p8-  /3p4+l 


dp  +  J  Jill  p'2  dp  +  O(p^) 

Pmax 


(5.2) 


Therefore  a  very  accurate  evaluation  of  I  can  be  obtained  by  using  the  trapezoidal 

rule  numerical  integration  for  the  first  term  in  the  right-hand  side  and  adding 

to  it  the  second  term  contribution  of  *^3/(2p  ).  A  step  size  Ap  -  0.02  was  used 

max 

for  the  numerical  integration  throughout.  Smaller  values  for  Ap  were  found  to 
change  the  numerical  results  only  in  the  seventh  or  eight  significant  figure. 

The  numerical  evaluation  of  p(0)  provided  agreement  to  eight  digit  accuracy 
with  the  analytical  value  given  in  section  III.  Furthermore,  p(0)  and  t(0)  were 
evaluated  using  both  of  the  expressions  valid  for  x  -  0  -  and  x  ■  0  +.  The  two 
results  coincided  up  to  the  seventh  digit,  again  showing  very  good  accuracy. 


The  results  for  x  and  p(x)  corresponding  to  a  wide  range  of  slot  lengths  x 
w  s 

are  presented  In  figures  2  and  3,  where  all  the  features  anticipated  In  Section 

III  can  be  easily  observed.  In  particular,  it  is  interesting  to  notice  the 

asymptotic  tendency  of  the  results  for  finite  xg  toward  that  of  the  limiting  case 

for  x^*.  Also  of  interest  is  the  fact  that  the  present  results  are  qualitatively 

very  similar  to  those  obtained  by  Smith  and  Stewartson  [4]  for  the  supersonic 

mainstream  case.  A  noticeable  difference  is  that  the  supersonic  mainstream 

-  3  x  <x‘xs) 

pressure  ahead  of  the  slot  [p  ■  (e  -  e  )]  is  essentially  independent 

of  the  slot  length  x  for  x  >5.  In  the  subsonic  mainstream  case  the  pressure 

s  s 

rise  ahead  of  the  slot  is  spread  out  over  a  much  longer  region  and  it  changes 
much  more  significantly  as  xg  grows  larger.  This  is  consistent  with  the  fact 
that  in  subsonic  maintream  flow  disturbances  usually  can  extend  upstream  over 
a  longer  region.  Incidentally,  note  that  the  graphical  results  in  Reference  4 
are  scaled  with  respect  to  the  p(0)  and  t(0)  values,  respectively,  so  that  p(x) 
and  t(x)  appear  to  be  independent  of  x  ,  which  is  untrue.  It  should  also  be 
pointed  out  that  very  similar  results  for  the  supersonic  and  subsonic  mainstream 
cases  were  found  by  Smith  [11]  for  the  problem  of  flow  past  a  small  hump  on 
a  flate  plate.  For  that  case  Napolitano  et.  al.  [8]  also  found  that  such  a 
similarity  was  maintained  by  the  nonlinear  results,  and  actually  for  both  cases 
separation  first  occurred  at  the  downstream  side  of  the  hump  as  anticipated  by 
linear  theory.  For  the  supersonic  slot  injection  problem  Napolitano  [6]  has 
found  that, due  to  nonlinear  effects, separation  first  occurs  ahead  of  the  slot 
for  a  length  of  the  same.x^  »  5.  An  even  more  remarkable  result, for  the  present 
subsonic  mainstream  case, is  discussed  next. 
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b)  Nonlinear  Equations  Results 

Solutions  to  the  nonlinear  equations  are  presented  for  one  slot  length 
only  xg  *  5.  A  comparison  can  then  be  made  with  the  supersonic  mainstream 
results  of  Napolitano  [6],  particularly  with  respect  to  the  occurrence  and 
location  of  separation.  A  small  value  of  the  wall  injection  velocity 
(V  *  0.1)  was  chosen  first,  in  order  to  check  accuracy  and  reliability 
of  the  numerical  technique  by  comparison  with  the  linear  theory  results. 

Both  the  numerical  technique  of  References  [8,  9]  with  the  appropriate 
boundary  conditions,  as  described  in  section  IV,  and  the  new  one  using  a 
variable  grid  and  a  direct  inversion  of  the  second  sweep  equation  were 
used  for  this  check.  The  wall  shear  results  obtained  by  the  two  different 
programs  are  shown  in  figure  4  together  with  the  linear  results.  The 
agreement  between  linear  and  nonlinear  solutions  is  quite  satisfactory, 
considering  that  the  nonlinear  solutions  are  only  first  order  accurate 
in  the  longitudinal  direction  x  and  use  a  fairly  large  step  size  Ax  »  0.25; 
also,  for  -  0.1  linear  theory  could  have  some  error.  Moreover,  the 
agreement  between  the  two  nonlinear  solutions  is  almost  perfect,  considering 
the  several  differences  between  the  two  schemes.  For  the  first  one,  for 
example,  the  step  size  Ax  was  constant  everywhere  (Ax  ■  0.25),  x^  and  x^ 
were  equal  to  -10  and  15  respectively,  and  the  second  sweep  equation  was 
modified  according  to  the  way  described  in  Reference  8,  section  IVc.  The 
new  scheme,  instead,  used  step  size  Ax  variable  outside  the  slot,  according 
to  the  stretching  provided  in  Reference  6.  x^  and  x^  were  equal  to  -20.8 
and  25.8  respectively  and  the  Cauchy  Integral  of  a  discontinuous  pressure 
gradient  is  evaluated  according  to  the  method  described  in  Reference  9, 
section  IV  3.  The  pressure  profiles  were  also  found  to  be  in  close 
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agreement.  A  step  size  study  was  performed  for  the  new  numerical  scheme. 

The  results  for  the  shear  and  pressure  at  the  leading  edge  of  the  slot 

are  shown  in  Figure  5;  the  second  order  accuracy  in  the  normal  direction  y 

and  the  first  order  accuracy  in  the  longitudinal  direction  x  are  clearly 

verified.  The  new  technique  was  then  employed  to  obtain  numerical  solutions 

for  larger  values  of  the  injection  velocity.  The  wall  shear  distributions 

obtained  for  three  values  of  V  ,  namely  0.5,  1  and  1.5,  are  presented  in 

figure  6.  These  results  are  very  interesting  and  quite  surprising. 

As  in  the  supersonic  mainstream  case  [6]  the  wall  shear  decreases  monotonically 

upstream  of  the  slot  and  has  a  dip  immediately  after  the  front  of  the  slot. 

It  then  decreases  again  along  most  of  the  slot  length.  In  the  supersonic 

free  stream  problem,  as  increased  to  a  value  of  about  1.1  the  minimum 

shear  became  negative  just  ahead  of  the  slot  indicating  the  occurrence 

of  separation  and  a  small  region  of  negative  shear  was  found  ahead  and 

through  the  slot  frong  for  Vw  *  1.5.  In  the  present  problem  the  shear 

dip  ahead  of  the  slot  also  increases  with  and  separation  would  appear 

likely  there  for  a  large  enough  V  value.  However,  in  this  case,  a  new 

w 

feature  appears,  i.e.,  there  is  a  more  rapid  shear  decrease  downstream  of 
the  slot  as  V  Increases.  This  is  small  and  very  localized  for  V  *  0.5 
but  it  becomes  more  and  more  significant  as  V '  increases.  For  V w  «  1.3, 
the  shear  actually  goes  to  zero  downstream  of  the  slot  and  for  ■  1.5 
a  reverse  flow  bubble  is  present  over  a  fairly  large  distance  on  the  plate. 

All  the  results  have  been  obtained  using  a  constant  step  size  in  the 
normal  direction  y.  Ay  ■  0,3  and  a  variable  one  in  the  longitudinal  direction 
x.  Ax  _>  0.25,  The  values  for  x^  and  x r  where  taken  to  be  -  20,8,  -  30,  -  35.2 
and  25.8,  35,  45.8  corresponding  to  Vw  <  0,5,  Vw  ■  1  and  Vw  -  1.5,  respectively. 


The  outer  edge  boundary  conditions  were  imposed  at  y  ■  12,  15,  18  and  21 
corresponding  to  Vw  -  0.1,  0.5,  1  and  1.5.  In  all  cases  the  solution  did 
not  change  appreciably  when  the  domain  of  integration  was  Increased.  Two 
points  are  of  special  interest.  First,  due  to  the  direct  inversion  of  the 
second  sweep  equation  the  optimal  time  step  for  convergence  was  found  to  be 
equal  to  two.  This  is  the  same  as  in  the  supersonic  mainstream  case  [6]. 

(For  the  technique  of  References  (8,9)  it  was  equal  to  one).  Second,  the 
convergence  rate  here  was  comparable  to  that  for  the  supersonic  mainstream. 

For  Vw  <  1  the  average  variation  for  6  over  one  complete  time  step  became  of 
the  order  of  10  7  in  less  than  one  hundred  iterations.  For  Vw  *  1.5  the 
convergence  rate  was  much  slower  so  that  it  had  to  be  reduced  to  one.  In  all, 
500  iterations  were  required  for  a  total  computational  time  of  about  8  CPU,  on 
the  AMDAHL/470  computer.  The  solution  was  initialized  with  the  uniform  shear 
flow  configuration  (u  -  y,  v  ■  6  ■  dp/dx  *  0  every  where) .  The  injection 
velocity  ^  was  started  at  the  value  of  0.1  and  increased  by  0.1  after  each 
time  step  until  the  final,  desired  value,  was  reached. 
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Figure  1.  THE  STRONG  SLOT  INJECTION  PROBLEM 
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Figure  2.  LINEAR  PROBLEM  WALL  SHEAR  PROFILES 


LINEAR  SOLUTION 


NONLINEAR  AND  LINEAR  NALL  SHEAR  RESULTS 


